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We develop a generic theory for the resonant inelastic light (Raman) scattering by a conduction 
band quantum plasma taking into account the presence of the filled valence band in doped semi- 
conductor nanostructures within a generalized resonant random phase approximation (RPA). Our 
generalized RPA theory explicitly incorporates the two-step resonance process where an electron 
from the filled valence band is first excited by the incident photon into the conduction band before 
an electron from the conduction band falls back into the valence band emitting the scattered photon. 
We show that when the incident photon energy is close to a resonance energy, i.e. the valence-to- 
conduction band gap of the semiconductor structure, the Raman scattering spectral weight at single 
particle excitation energies may be substantially enhanced even for long wavelength excitations, and 
may become comparable to the spectral weight of collective charge density excitations (plasmon). 
Away from resonance, i.e. when the incident photon energy is different from the band gap energy, 
plasmons dominate the Raman scattering spectrum. We find no qualitative difference in the res- 
onance effects on the Raman scattering spectra among systems of different dimensionalities (one, 
two and three) within RPA. This is explained by the decoherence effect of the resonant interband 
transition on the collective motion of conduction band electrons. Our theoretical calculations agree 
well (qualitatively and semi-quantitatively) with the available experimental results, in contrast to 
the standard nonresonant RPA theory which predicts vanishing long wavelength Raman spectral 
weight for single particle excitations. 



PACS numbers: 73.20.Mf; 78.30.Fs; 71.45.-d 
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I. INTRODUCTION 



In recent years the elementary electronic excitation spectra of a variety of doped semiconductor nanostructures, such 
as two dimensional (2D) quantum wells (QW) heterostructures, superlattices, and more recently, one dimensional (ID) 
quantum wire (QWR) systems, have been studied extensively both experimentally [p]-pl| and theoretically [ p^ 
Rich experimental spectra of the elementary electronic excitations (such as charge density excitations (CDIC) 
density excitations (SDE), and single particle excitations (SPE), for both intrasubband and intersubband excitations) 
in these systems are typically experimentally investigated by using the resonant Raman scattering (RRS) technique, 
which is a powerful and versatile spectroscopic tool to study interacting electron systems. In the RRS experiment, 
external photons are absorbed at one frequency and momentum, uji and k^, and emitted at another, ojf and k/, 
creating one or more particle-hole pair (or collective) excitations in the conduction band. The energy and momentum 
difference between the incident photon and the scattered photon is Stokes shift, indicating the dispersion of the relevant 
elementary electronic excitation created in the system. In the so-called polarized RRS geometry with the incident and 
scattered photons having the same polarization, the excited electrons have no spin flips during the scattering process, 
which therefore corresponds to the elementary charge density excitations of the system. At low temperatures (which 
is of interest to us in this paper) there is no real absorption of elementary excitations by the incident photon and the 
anti-Stokes line is not of any importance. (We use % = 1 throughout this paper.) 

In the standard theory |p6|-p0[ , which ignores the role of the valence band and simplistically assumes the external 
photons to be interacting entirely with conduction band electrons, the polarized RRS intensity is proportional to 
the dynamical structure factor of the conduction band electrons and therefore has strong spectral peaks at 

the collective mode frequencies at the wavevectors defined by the experimental geometry. The dynamical structure 
factor peaks correspond to the poles of the reducible density response function, which are given by the collective 
CDEs (plasmons) of the electron system in the long wavelength limit. In particular, the single particle electron- 
hole excitations, which are at the poles of the corresponding irreducible response function, carry no long wavelength 
spectral weight (about three orders of magnitude weaker than the CDE spectral weight at the typical wavevector, 
10^ cm~^, accessible in RRS experiments) in the density response function (according to the /-sum rule |^^). The 
SPE therefore should noi, as a matter of principle, show up in the polarized RRS spectra in any dimensions. The 
remarkable experimental fact is that there is always a relatively weak (but quite distinct) SPE peak in the observed 
polarized RRS spectra in addition to the expected CDE peak. This experimental presence of SPE peak in RRS cannot 
be explained by the standard theory, which, however, does give the correct mode dispersion energy for both the CDE 
and the SPE, but fails to explain why the SPE spectral weight is strongly enhanced in the RRS experiments. This 
puzzling feature |p],p|Jl9|] of an ubiquitous anomalous SPE peak in addition to the expected CDE peak (or equivalently 
a two-peak structure) occurs in one, two, and even in three dimensional doped semiconductor nanostructures It 
exists in low dimensional semiconductor systems both for intrasubband and intersubband excitations. 

Many theoretical proposals [p^|- p"5|j2^j2^ have been made to explain this two-peak RRS puzzle. Ad hoc proposals 
P,p2[ have been made in the literature that perhaps a serious breakdown of momentum or wavevector conservation 
(arising, for example, from scattering by random impurities) is responsible for somehow transferring spectral weight 
from large to small wavevectors, because the usual linear response theory predicts that at very large wavevectors (an 
order of magnitude larger than the experimentally used RRS wavevectors) , where the CDE mode is severely Landau 
damped, the dynamical structure factor should contain finite SPE spectral weight corresponding to high energy 
electron-hole excitations. Apart from being completely ad hoc, this proposal also suffers from any lack of empirical 
evidence in its support — in particular, the observed anomalous SPE peak in the RRS spectra does not correlate with 
the strength of the impurity scattering in the system. We have recently systematically analyzed jl^] all the proposed 
mechanisms within the non-resonant RRS theory (i.e. without incorporating the valence band in the theory, assuming 
simply the inelastic light scattering process to be entirely confined to the conduction band free carrier system) leading 
to the conclusion that none of the proposed nonresonant mechanisms can explain the ubiquitous two-peak (the lower 
energy SPE peak and the higher energy CDE peak) structure of the observed RRS spectra. 

We have recently reported in a short letter a new resonant RRS theory by generalizing the nonresonant RPA 
theory to include the filled valence band in the semiconductor, reflecting the two-step resonant nature of the RRS 
process. The purpose of the current article is to provide the details of our resonant RRS theory, and more importantly, 
to present RRS results for 2D and 3D systems which have not been discussed earlier in the literature (our earlier 
letter |12 presented only ID RRS results). The observed experimental RRS phenomenology in ID, 2D, and 3D 
systems being very similar qualitatively, our generic interband-resonant RRS theory, as reported herein, provides the 
conceptual theoretical foundation for understanding RRS spectroscopy in doped semiconductor structures. 

In this context, we emphasize that the striking similarity of the experimental RRS spectra in one, two, and three 
dimensional semiconductor systems suggests that the problem (namely, the two peak nature of the RRS spectra with 
the conspicuous presence of the "forbidden" SPE peak) is not specific to ID systems, where our earlier theory llj] 
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was applied. The ubiquitousness of the strong SPE spectral weight in the RRS experiment (independent of system 
dimensionality, dependent only on the resonant nature of the experiment) suggests that the theoretical explanation for 
this puzzle must arise from some generic physics underlying RRS itself, and cannot be explained by the non-generic 
and manifestly system-specific theories which have been made occasionally in the literature. The resonant RRS theory 
presented herein (and in our short letter) provides a generic explanation for the two-peak structure of the RRS spectra 
by establishing that the so-called low energy anomalous SPE feature in the RRS spectrum arises entirely from the 
resonant two-step nature of the RRS experiment, and cannot be the explained within any non- resonant theory. 

In this paper we provide (within the resonant RPA scheme) the compellingly generic theory for RRS experiments by 
including the valence band electrons during the scattering processes for one, two and three dimensional semiconductor 
systems, following our earlier short paper |l2) on ID systems. We find that the RRS spectral weight at SPE energy is 
a strong function of the resonance condition — the SPE spectral weight is substantially enhanced when the incident 
photon frequency is near the semiconductor band gap resonance energy, and decreases drastically away from the 
resonance. It is important to emphasize that this feature of our theory agrees with experimental observations — the 
anomalous SPE peak exists only around resonance and its spectral strength decreases off resonance. Our results show 
similar qualitative behavior for the RRS spectra in one, two and three dimensional systems. 

One dimensional systems [ p4[ actually pose a special (and subtle) problem with respect to understanding the two- 
peak RRS spectra because ID electron systems are generically ||2^,Q Luttinger liquids (i.e. non- Fermi liquids) which 
have no quasiparticle (SPE) excitations whatsoever. The elementary excitations in ID electron systems are bosonic 
spinon and holon collective modes. It is therefore conceptually problematic to comprehend how an anomalous "SPE" 
feature can arise in ID semiconductor quantum wire RRS spectra as has been experimentally observed |l- 
The issue of understanding ID RRS spectra from a Luttinger liquid viewpoint has been recently discussed 
the literature, and we refrain from further discussing this point in this article since this is beyond the scope of our 
work. In particular, our use of generalized RPA enables us to develop a unified consistent theory of resonant RRS in 
arbitrary dimensions (including ID), and the Luttinger liquid nature of ID quantum wires is not of any relevance in 
our theory. We mention, however, that a complete Luttinger liquid theory of ID RRS experiments has recently been 
developed and this Luttinger liquid theory builds on the resonant nature of our work presented in this article. 

The rest of this paper is organized as follows: in Sec. II we describe the theory of nonresonant and resonant Raman 
scattering process in RPA; in Sec. Ill we present and discuss our calculated RRS results for one, two, and three 
dimensional GaAs semiconductor systems; we then summarize our work in Sec. IV. All the results shown in this 
paper are for GaAs-based systems, but obviously the theory applies to any direct band gap semiconductor material. 




II. THEORY 



In Fig. 1(a) we depict the schematic diagram | p^ , p7| -|2g[| for the two step process (steps 1 and 2 in the figure) 
involved in the polarized resonant Raman scattering spectroscopy at the Eq + Aq direct gap of electron doped GaAs 
system where an electron in the valence band is excited by the incident photon into an excited (i.e. above the 
Fermi level) conduction band state, leaving a valence band hole behind (step 1); then an electron from inside the 
conduction band Fermi surface recombines with the hole in the valence band (step 2), emitting an outgoing photon 
with an energy and momentum (Stokes) shift. The RRS process is a two-step process involving steps 1 and 2 with 
the net result of there being an elementary electron excitation created in the conduction band through intermediate 
valence band states as shown in Fig. 1. The non-resonant approximation to RRS ignores the intermediate valence 
band states and approximates the RRS process to be taking place entirely within the conduction band of the system, 
as shown by the step 3 in Fig. 1. The whole point of the theory developed in this paper is that the non-resonant 
step 3 is not equivalent to the resonant scattering involving steps 1 and 2. Note that the resonant process involving 
steps 1 and 2 explicitly depends on the incident photon energy, while the non-resonant approximation depicted in step 
3 depends only on the energy difference between the incident and the scattered photons and not on the incident photon 
energy. This difference turns out to be crucial in the RRS theory, and the resonance condition in the incident photon 
energy gives rise to the anomalous SPE-like feature in the RRS spectra as shown below. Electron spin is conserved 
throughout the scattering processes since we are considering only the polarized geometry where no spin flip occurs. 
As mentioned before we use the random phase approximation (RPA) in our theory taking care to generalize it to the 
resonant situation involving steps 1 and 2. In the RPA one neglects all exchange-correlation effects (e.g. self-energy 
and vertex corrections due to electron-electron interaction), including only the long range Coulomb interaction Vc(q) 
in the dynamical screening by the electron system so as to correct the noninteracting irreducible response function 
to the reducible response function. Following a preliminary discussion of the Coulomb interaction in 1, 2 and 3 
dimensional semiconductor system in subsection A below, we then develop the nonresonant and the resonant RRS 
theories in sections B and C respectively. Our theory is entirely within the effective mass approximation, and we 
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parameterize the electron system in the semiconductor by electron (me) and hole (m/J effective masses corresponding 
to the top (bottom) of the conduction (valence) band and by a background lattice dielectric constant Eq. 



A. Coulomb Interaction 



The realistic (bare) Coulomb interaction in the artificially confined semiconductor nanostructures depends strongly 
on the confinement geometry of the systems. In the bulk 3D semiconductor materials, the unscreened Coulomb 
interaction has a long-ranged 1/r decay in the real space, and has the following Fourier transform in momentum 
space: 



' So J |r| £o|qP' ^ ' 



where e is the electron charge and Eq is the dielectric constant of the background material (about 12 in the GaAs 
semiconductor system). We use the static (eg) lattice dielectric constant in our theory rather than the more conven- 
tional high frequency (eoo) dielectric constant in defining the Coulomb interaction, K:(q), because inclusion of eq is 
known to approximately account for the polaronic electron-phonon interaction in the system, which is rather weak in 
GaAs because of its low Frohlich coupling constant 0.07). In the 2D semiconductor quantum well system, modern 
fabrication techniques have produced very narrow 2D wells (of nanostructure size < 100 A in GaAs in the confinement 
direction), leading to an almost pure 2D electron system. It is therefore a good approximation to assume the well 
width to be zero in our calculation, giving a 2D Fourier transform of the Coulomb interaction: 

K^°(q) = - / ??e-- = ^. (2) 



£o J \r\ Eo\q\ 

Inclusion of the confinement wavefunction effect in the theory is straightforward and leads to a form factor /(q) 
(< 1) multiplying l/?^(q) in the theory. For the ID semiconductor quantum wire system, we have to consider the 
realistic finite width of the wire (i.e. the relevant ID form factor effect) because the ID Fourier transform oi 1/r 
potential (i.e. / fire"^'"/|r|) diverges logarithmically requiring regularization by a length cut-off associated with the 
typical confinement size. Therefore the Coulomb interaction for the finite width quantum wire is obtained by taking 
the expectation value of the 2D Coulomb interaction (assuming the width in the z direction to be zero for simplicity 
as in our 2D model in Eq. 2) over the confinement wavefunction along the transverse direction (y) of the wire. We 



then have |30| the following Coulomb interaction matrix element in the ID QWR structure of finite width: 



2 foo j-ca 

'-■0 J —OO J —OQ 



= — dydy'\My)\'\<l>jiy')\'Ko{q\y~y'n, (3) 

^0 J-oo 

for interaction between electrons of subband i and subband j. 4>i(y) is the electron wavefunction of ith subband of 
the QWR along the transverse direction. In this paper we assume that only the lowest {i = 1) ground conduction 
subband is occupied (i.e. subband spacing AE12 > Ep at zero temperature and all the higher energy subbands 
are empty) and neglect any intersubband transition, so that the subband index i = 1 throughout and will not be 
explicitly shown. Kq{x) in Eq. (3) is the zeroth-order modified Bessel function of the second kind. The exact form 
of wavefunction (/)(j/) depends on the confinement geometry of the QWR system. For simplicity we assume the QWR 
confinement potential to be the ID infinite square well in the y direction. This turns out to be a good approximation 
for the electrostatic gate-controlled confinement in the presence of the self-consistent Hartree potential due to the free 
electrons themselves BG]. The the confinement wavefunction 4'{y) is (for the ground subband): 




if-a/2 < y < a/2 
otherwise. 



where a is the wire width in the y direction. Using Eqs. (3) and (4) we can numerically calculate the effective ID 
Coulomb interaction for the semiconductor QWR system. Unlike the power-law behavior of Coulomb interaction 
in the higher dimensions (Eqs. (1) and (2)), V™{q) has a weak logarithmic divergence, —2e'^h\{qa)/EQ, in the long 
wavelength limit [q ^ 0). Because of this logarithmic dependence of Vc{q) on q (as 5 ^ 0), the precise value of the 
wire width (a) is not particularly important in our theory, making our simple infinite square well approximation a 
reasonable one for our purpose. 
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B. Nonresonant Raman scattering 



In the presence of an external photon field the interacting Hamiltonian between the free electron gas and the 
radiation field is assumed to be obtainable from the standard gauge-invariant prescription p ^ p — eA/c, 

where A is the radiation field (photon) vector potential operator and c the speed of light. The Hamiltonian including 
the radiation field and the electrons (i.e. the free carriers induced by doping) in the semiconductor conduction band 
can therefore be written as (we neglect the spin-photon interaction considering only polarized RRS spectra where 
spins do not play any explicit role) 



H = (p« - eA(x,, t)/cf + V V;(x, - Xj) 



N 

E 



2 ^ 



Ho 



N 



p, • A(x,, t) + 2 

nieC 2,m,,c 



Hi 



A(x,,i) 



(5) 



in the effective mass approximation with me being the effective electron mass of the semiconductor conduction band. 
We have made the transverse gauge choice V • A(xi, t) = 0, for the radiation field, leading to Pi ■ A = A ■ pi as 
used in Eq. (5). is the Hamiltonian of electrons interacting with Coulomb potential without the radiation field, 
and Hi is the electron-photon interaction Hamiltonian which plays a crucial role in the Raman scattering problem. 
Figs. 1(b) and 1(c) correspond to the scattering processes induced by the linear (p • A) term and the quadratic (A^) 
term respectively in the second quantization representation. The p • A term creates and annihilates one photon in the 
state it acts on, having no contribution to the scattering rate in the first order time-dependent perturbation theory 
since there is no net change of photon numbers. The quadratic A? term, on the other hand, gives a non-vanishing 
first order contribution to the scattering rate because photons are created and annihilated at the same time in such 
scattering processes as shown in Fig. 1(c). In principle the second order contribution of p • A term in the time- 
dependent perturbation theory is of the the same order as the first order contribution from the A^ term as a simple 
power counting in the coupling constant e/c shows. This second order contribution, which plays a role in the RRS 
phenomenon, will be studied and discussed in more details in the next section. We can simply neglect this (p • A) term 
in Hi if we are interested only in the nonresonant Raman scattering regime, either because the incident photon energy 
uji is off-resonance i.e. far from the direct band gap, (~ 1.5 eV in GaAs), or because we only want to consider a 
nonresonant process as in step 3 in Fig. 1(a). The A? term, being a scalar field operator which commutes with the 
electron field, V'(x), leading to the perturbative Hamiltonian, Hi (neglecting the p • A term) being proportional to 
the electron density operator, n(x) = ■0j(x)^/'s(x). The nonresonant (corresponding to the step 3 process in Fig. 
1) Raman scattering intensity at frequency shift uj and momentum transfer q therefore can be calculated from the 
dynamical structure factor (the imaginary part of the density response function) in the linear response theory pO|,pl[ : 



dQdio 



-ImH(q, oj) = Im 



dte^'^'{[n\q,t),n{q,0)])o 



(6) 



where 



)o is the ground state expectation value, and n(q, t) is the electron density operator, n(q, t) = 
{t)cy:,s{t), with c]^^(ck,s) the electron creation(annihilation) operator for momentum k and spin s. In 
the standard many- body theory, this (reducible) response function can be obtained by the reducible set of polar- 
ization diagrams l2^,^ltl (Dyson's equation, see Fig. 2) formed by the irreducible conduction band polarizability, 
Hq (q, uj) , for the scattering process where one has an electron and a hole in the conduction band: 



H(q, u) = no(q, lo) + Ho(q, lo) y,(q)Ho(q, u) 
^ no(q,cj) _ Ho(q,^) 
" l-V;(q)Ho(q,c^) " e(q,c^) ' 



(7) 



where e(q, w) = 1 — Vc{q)Tlo{(l,uj) is the dynamical dielectric function. 

In the random phase approximation (used in this paper), the irreducible polarizability, Ho(q, w), is approximated 
by the noninteracting electron- hole bubble, H|j'^^^(q, a;), without any self-energy or vertex correction. RPA is known 
to be a good approximation Jl6|- |T8| , p7| -p9| in two- and three-dimensional electron systems for calculating plasmon 
(or CDE) properties. It is also a good approximation for collective mode dispersion in one-dimensional electron 
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systems and gives a ID plasmon dispersion which agree with the exact Luttinger Uquid theory |Q. The expression 
of nj^^'^(q, w) for a d-dimensional system is 



—2i f 

n^P^(q,cj) = j^^^ J dvdp Go(p,i^)Go(p + q,i^ + w) 



{2ttY J UJ + - p2/2me + (p - q)2/2me 

where Go(p, ly) is the bare conduction band electron Green's function and no(p) = 0{kF — |p|) is the zero temperature 
noninteracting momentum distribution function of conduction band electrons. 7 is a phenomenological damping term 
associated with impurity scattering (or other broadening mechanism), which is taken to be small (7 ^ Ep) in our 
numerical calculation. The damping term, 7, introduces finite widths to the spectral peaks in the dynamical structure 
factor of Eq. (6), but does not affect the peak position and spectral weight in any significant method. The imaginary 
part of the irreducible polarizability IlQ{q,uj) (which is now approximated by n|^^^(q, a;) in our paper) gives rise to 
the single particle excitation, which is typically very small at long wavelengths due to the dynamical screening effect 
of Eq. (7). In Fig. 3 we show as shaded regions the SPE continua (where Imno(q, w) 7^ 0) within RPA for one, two, 
and three dimensional systems. Note that, in contrast to 2D and 3D systems, the ID SPE continuum is very restricted 
in the long wavelength limit (q <C fc^?). In higher dimensions, the SPE continuum is gapless for any finite wavevector 
smaller than 2fci?, but it is gapped in ID due to energy-momentum conservation induced phase space restriction. 
Using Eqs. (6)-(8) we can calculate the nonresonant Raman scattering spectra and the plasmon (CDE) dispersion 
(shown in Fig. 3) to compare with the experimental results and the resonant theory results discussed below. The 
calculated spectra are shown in Figs. 4(a)-(c) for one, two and three dimensional systems respectively. We discuss 
these results in details in Sec. III. 



C. Resonant Raman scattering 

We now consider the full resonance situation (step 1 and 2 in Fig. 1) including the valence band which obviously 
],^ -14 3^ plays a crucial role in the RRS experiment because the external photon energy must be approximately 



equal the Eq + Aq direct gap for the experiment to succeed. In the RRS process the incident photon is absorbed and 
a scattered photon with the appropriately shifted frequency (and wavevector) is emitted. Electron spin is conserved 
throughout the scattering process. As discussed above, there are two steps (steps 1 and 2 in Fig. 1(a)) involved in 
the polarized RRS spectroscopy and both of these two steps of inelastic scattering result from the p • A term of Hj in 
Eq. (5) (see Fig. 1(c)). When the incident photon frequency is equal to the direct band gap energy, Eq, the second 
order "resonant" perturbative contribution of the p • A term becomes important and comparable to the first order 
contribution of A^ term, leading to an electron interband transition between the conduction band and the valence 
band. The interaction Hamiltonian of the RRS theory with external photon momentum, k, and frequency, w, can be 
expressed in the second quantization representation as 

p,cr 

+ e^"* E [4- (0«P+k,- (t) (Ocp+k,. it)] , (9) 

p,cr 

with Ck,cr and Wk,(T being the annihilation operators of conduction band and valence band electrons respectively. The 
electron-photon coupling vertex, p ■ e (where e is the light polarization), has been assumed to be constant for 
simplicity. Applying the time-dependent perturbation theory to the ground state, |0), characterized by a conduction 
band Fermi sea and no holes in the valence band (at zero temperature), we have the following transition amplitude 
from ground state, |0), to the nth excited state, \n): 

c„iT)^r\hf' di2(n|ff/'-'(ii)i?^'"'(t2)|0), (10) 

J-T/2 J-T/2 

where we have changed the time integration range from the conventional {0,T} to {— T'/2,T/2} for the convenience 
of changing variables later. By substituting the explicit form of Hj and choosing the specific channel of backward 
scattering (the so-called back-scattering geometry), kj = — k/ = q/2 and ujij — U±uj/2, without any loss of generality, 
we obtain the transition rate (ignoring excitonic and self-energy effects), W, to be 
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W= lim - 



= lim - V / dh / di2 / dt[ / rf4ein(*2-«'i+ti-*2)gia,(ti+ti-ti-t.)/2 

T^oo T J -Til J -Til J-T/2 J-T/2 



Pl ■ ■P4 
(71 ■■■(T4 



T/2 .ii .T/2 .t'l 

dtl / dt2 / rft'i / 
T/2 J-T/2 J-T/2 J-T/2 

><Hi-q/2,<Ti(*2)^^P2,<r2(i'l)43,a3(*l)^P4-q/2,a4(*2))0 

x(Cpi,ai(4)Cp,_q/2,<.,(*'l)Cp3-q/2,a3(tl)c^4,a4(*2))0. (H) 

Since the valence band is completely filled in the ground state at zero temperature, we have only one contraction of 
the valence band electron operators, which is assumed to be noninteracting for simplicity: 

(^Pl-q/2,<Tl(*2)^^P2,-T2(*'l)w],3,<,3(tl)«p4-q/2,a4(*2))0 

= (^ii-q/2,<Ti(4)^^P2,<T2(ii))o(wJ,3,,r3(^l)^'p4-q/2,a4(*2))o. 

- S ^ ,nS S ,nS e'^vipi-l/2){t'i-t'2) iEUp3)itl-t2) (,2) 

~ ''pl,P2+q/2'Vi,CT2'-'p3,P4-q/2"<T3,CT4'^ ^ ) \^'^) 

where Ey{p) = — p^/2m^ is the kinetic energy of the valence band electrons. Setting new implicit time variables 
(ti,2 ii± ti/2 and t'l 2 — > t2 ± ^2/2) and using the quasi-particle approximation for the electron operator, Cp,cr(i — 
t/2) = Cp,o-(f)e~*^°(P^*/^, we can obtain the transition rate, after evaluating the ti and t2 integrals, 

1 _ rT/2 _ 

W= lim - / dh / d<2e*'^(*=-*i) 

T->oo T J_T/2 J-T/2 

X E A*(pi,q)^(p2,q)(Cp,+q/2,<.i(t2)cJ^_q/2,ai(^)Cp2-q/2,a2(^'l)ci,+q/2,<.,(^'l)>0 

Pl ,P2 
(7l ,0-2 

/•OO 

= / dte'-*(iVt(q,f)Ar(q,0))o, (13) 
where the resonant "density" operator, iV(q, t), is defined to be 

t) = -$Z^(P'q)Cp-q/2.-(^)4+q/2,-(^) 

p,(T 

= E A(p,q)c;^q/2,.(*)cp-q/2,.(0 (14) 

p,cr 

for q 7^ with the matrix element A(p, q): 



^(p, q) = ^ 



Eg-n + (E^ip - q/2) + Ec{p + q/2))/2 - E^p) + iX 

1/Ef 



Eo. + {l+ 0(P' - 1) + qV4 + i>^/EF ' 



(15) 



Here E^ = Ef~^ {Eg + {1+£_)Ef — with ^ = rric/my: p = p/kF] q = q/fc^, and Ep = Ec{kp) is the Fermi energy of 
the conduction band electrons. A is a phenomenological broadening factor we introduce to include roughly all possible 
broadening effects, e.g. finite imaginary part of electron self-energy (the quasi-particle life time), finite impurity or 
disorder scattering, and any broadening or damping arising intrinsically from the photon field or the associated optical 
scattering. We take A to be small {=0.02Ef) in the numerical calculation. Note that the phenomenological parameter, 
A, is a resonance broadening parameter (associated with the band to band process) to be contrasted with the simple 
spectral broadening parameter, 7, of Eq. (8), which is purely a conduction band phenomenological parameter. In 
our leading order RRS theory A (of Eq. (15)) and 7 (of Eq. (8)) are completely independent phenomenological 
relaxation or damping terms (both of which should be small, 7 and A <C Ep, for our leading order theory to be 
sensible). Calculation of 7 and A is beyond the scope of the leading order theory — it is entirely possible that in a 
more complete theory including quasiparticle self-energy and vertex corrections as well as electron-impurity scattering 
and the electron-photon interaction, 7 and A will tern out to be related. 
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Comparing Eqs. (13)-(14) with Eq. (6), we find that the effect of resonance (i.e. photon induced interband 
transition) on the conduction band electrons is the matrix element A{p,q), which arises from the time difference 
between the excitation of one electron from valence band to the conduction band (step 1) and the recombination 
of another electron from inside the conduction band Fermi surface with the hole in the valence band (step 2). The 
resonance condition is parameterized by the dimensionless parameter, i?^, with Ei^ = being the precise resonance 
condition. In the following discussion we define "off resonance" as \Ei^\ ^ 1 and "near resonance" as \E^\ <C 1. Off 
resonance the spectral weight decreases as as can be seen from Eq. (15). Near resonance the singular properties 

of ^(p, q) enhances the spectral weight nontrivially. The calculation of the RRS spectrum is therefore reduced to 
the evaluation of the correlation function of Eq. (14), which in the resonant RPA approximation (i.e. neglecting all 
vertex correction of the irreducible polarizabilities, see Fig. 5) is obtained to be 



W « -Im 



nr^(q,c.) + 



nfPA(q,c.)nr^(q,c.)F,(q) 



(16) 



where 



n^^^(q,c.) = 



dp 



lAjp, q)P(no(p + q/2) - no(p - q/2)) 
LU + i-f-Ecip + q/2) + Ecip - q/2) ' 



(17) 



and 



ttRPA/„ ^_ -2 f ^(p,q)(no(p + q/2)-no(p-q/2)) 



uj + i-f-Ec{p + q/2) + Ec{p- q/2) ' 



(18) 



f,nPA,,._ -2 f , >l*(p,q)(no(p + q/2)-no(p-q/2)) 

The dynamical dielectric function, e(q, oj), is the same as defined in Eq. (7) within the same RPA formulae (Eq. (8)). 
Note that resonance effects arising from A{p, q) (i.e. considering the full two step process involving both conduction 
and valence bands rather than just the effective single step process (step 3 of Fig. 1(a)) within the conduction band) 
are nonperturbative and depend crucially on the exact value of the incident photon energy. In the nonresonant theory, 
by contrast, the incident photon energy does not enter into the calculation of the spectra, only the frequency shift lo 
matters. 



III. RESULTS AND DISCUSSIONS 



In Figs. 3 and 4 we show the energy dispersion and the dynamical structure factor respectively of the nonresonant 
RRS spectra in the RPA theory for ID, 2D, and 3D semiconductor GaAs systems. We emphasize that all earlier 
theoretical works on RRS spectroscopy, with the only exception of our earlier brief communication | [T^ , use the 
nonresonant approximation. The sold lines in Fig. 4 are the RRS spectrum profiles in the long wavelength limit (small 
momentum transfer, \q\—0.1kp), while the dashed lines are the results of larger momentum transfer for comparison. 
(The experimental situations correspond to the long wavelength limit, with |q| <^C kp.) Two elementary excitations are 
observed in the nonresonant spectra (Fig. 4) at two separate peaks: one is the single particle excitation at lower energy 
and the other one is collective charge density excitation at higher energy. (Note that we use very small damping, 
7 = 10~^ Ep, in Fig. 4 in order to resolve the small SPE weights; larger 7 smears out the SPE continuum completely.) 
We first mention that the RPA calculated energy dispersions of both modes (SPE and CDE) agree quantitatively 
with the experimental RRS results |ll|, |l6| , p^l 9y24|j3^ . However, the theoretically calculated nonresonant dynamical 
structure factor in Fig. 4 is entirely dominated by the collective CDE mode, the SPE mode, while being present in the 
results, carries negligible and unobservable spectral weight. This is entirely inconsistent with the "two-peak" structure 
observed in the experimental RRS spectra [|l| where the two peaks carry comparable spectral weights. In the large 
momentum transfer results (which are outside the experimentally accessible regime) shown in Fig. 4 (dashed lines), 
one finds that SPE spectral weights are somewhat enhanced over the long wavelength results, and correspondingly 
CDE weights decrease for large momentum scattering due to the strong Landau damping of plasmons (CDE) to single 
particle excitations which become allowed at large wavevectors. The SPE spectral weight is still much weaker (by three 
orders of magnitude) than the CDE weight even at large wavevectors, and in addition, the incoherent SPE continuum 
is severely broadened in this large momentum scattering channel. Note that this situation (i.e. negligible theoretical 
spectral weight at SPE) does not change ||jl^,^,Q even if one goes beyond RPA and includes vertex corrections 
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(e.g. Hubbard approximation or time-dependent local density approximation) in the irreducible response function. 
Therefore, as long as resonance effects are neglected (and thus one includes only the step 3 of Fig. 1(a) ignoring 
the interband resonance process), the calculated RRS spectra at experimentally accessible wavevectors produce only 
observable CDE peaks in contrast to the experimental two-peak situation which, in addition, always finds at resonance 
SPE spectral weight to be comparable to the CDE spectral weights |l|- [Tl|JT9| , p^ . The nonresonant theory is therefore 
in qualitative disagreement with experiments as it fails to account the observed two-peak RRS spectra. 

In Figs. 6-7 we show our results for the polarized RRS spectroscopy of the same ID, 2D and 3D systems as in Fig. 4 
within the resonant RPA theory (Eqs. (13)-(19)) in the long wavelength region (|q|=0.1 kp). RRS spectra for different 
resonance conditions, i.e. for different values of Ei^, are shown in Fig. 6 with a larger value of the impurity broadening 
parameter (7 = O.ObEp, fifty times greater than the 7 used in Fig. 4) in order to compare with the experimental 
RRS profiles. The lower(higher) energy peak is associated with the SPE(CDE) of the electron systems. The most 
important qualitative feature of the resonant theory results is the great enhancement of the SPE spectral weight 
compared with the nonresonant theory. The three figures. Figs. 6(a), 6(b), and 6(c) (corresponding to the results of 
ID , 2D, and 3D systems respectively) have qualitatively very similar behaviors: (i) the overall spectral weights decay 
very fast off resonance (i.e. for large |i?ij|); (ii) the peak positions of the SPE and CDE in Fig. 6 are the same as 
the nonresonant excitation energies in Fig. 4, i.e. resonance does not affect the energy dispersion of the elementary 
electronic excitations; (iii) the spectral weight of SPE (lower energy peak) is essentially zero far away from resonance 
{\Euj\ > 0.2) where the CDE (higher energy peak) dominates similar to the nonresonant spectra in Fig. 4 (except for 
the larger value of 7 used in Fig. 6); (iv) near resonance {\Ei^\ < 0.2) the SPE spectral weight is greatly enhanced — in 
fact, the SPE spectral weight becomes comparable to or even larger than the CDE spectral weight, in sharp contrast 
to the nonresonant theory (where the SPE weight is always extremely small at long wavelength). In Fig. 7 we plot 
our calculated RRS spectral weight ratio of CDE/SPE as a function of the resonance condition, explicitly showing 
the dramatic effect of resonance on the SPE spectral weight. We emphasize that this spectacular enhancement of 
SPE spectral weight in the full two step resonant scattering process (over the simple one step nonresonant effective 
theory) is a nonperturbative effect in our theory. Our calculated spectra at resonance are in excellent qualitative 
agreement with the corresponding experimental RRS spectra shown in Ref. , where the SPE spectral weight dies 
off rather quickly as the incident photon energy goes off-resonance. From our results presented in Fig. 7, we also find 
that the spectral weight ratio of the CDE to the SPE has very similar resonance behaviors for systems of different 
dimensionalities, consistent with the experimental findings and indirectly ensuring the validity of RPA theory in the 
RRS spectroscopy, at least in the experimental parameter regimes. 

To understand the resonance condition dependence (on E^^) of Fig. 7, we should explain the resonance effects 
not only on the SPE continuum, but also on the CDE modes around the resonance region. In some sense the ex- 
treme resonance condition, Ei^ = 0, may be thought of as providing an indirect mechanism for the breakdown of 
the wavevector conservation for the scattering process considered only within the conduction band in the prevailing 
nonresonant theory where the virtual valence band effects are ignored (i.e. step 3 in Fig. 1(a)) — thus our theory 
preserves the essence of the "massive" wavevector breakdown mechanism proposed in Ref. but in a very indirect 
sense because no impurity scattering is involved. Instead, participation by the valence band introduces the effective 
mechanism for wavevector conservation "breakdown" through virtual interband process not included in the nonreso- 
nant theory. In particular, the function j4(p, q) defined in Eq. (15) provides the "wavevector conservation breaking" 
mechanism by mixing conduction and valence band wavevectors non-trivially; if A(p,q) is a constant, there is no 
resonant enhancement of the SPE mode. Equivalently, the dependence of A{p, q) on two different wavevectors is 
the effective wavevector conservation breakdown mechanism. Mathematically we can start from the RPA dynam- 
ical structure factor defined in Eq. (16), where the CDE spectral weight is given by the numerator of the second 
term, ^^(q, w)nf ^^(q, w) Fc(q), at the CDE dispersion energy determined by the zero of the dielectric function 
(e(q, cj) = 0). Off-resonance, the function A(p, q) is just a slowly varying function of momentum, p, in the integral 
range \p ± kp\ < q/2 obtained by the occupancy factor, no(p -I- q/2) — no{p — q/2), in Eq. (17)-(19), and therefore 
the RRS spectra show behavior (i.e. the CDE dominance the SPE) similar to the standard RPA results (see Fig. 
4) except for the overall decreasing weight factor, E~'^. Near resonance {E^^ ~ 0), however, the resonance function 
A{p,q) in Eqs. (18) and (19) can essentially cancel the contribution from the other integrand in the polarizabilitics 
(due to its sign change at |p| = hp), so that the CDE spectral weight (coming essentially from the second term in 
Eq. (16)) cannot be as strongly enhanced by resonance as the SPE weight, which arises mostly from the Il2^^{q,uj) 
in Eq. (17). Therefore the sign change of the resonant function, A{p,q), is responsible for the relatively weaker 
enhancement of the CDE weight compared to the SPE weight near resonance. We note that Eq. (16), defining the 
resonance spectral weight in our theory, has two terms, both of which are important in giving rise to a strong SPE 
spectral feature in the RRS spectra under resonance conditions. 

Finally we give a simple explanation for the breakdown of Luttinger liquid theory in the ID RRS process near 
resonance. It is well known that ID electron systems are best understood as a Luttinger liquid, where collective 
excitations are the only possible excitations and no single particle excitations exist, for the conduction band electrons. 
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However, Luttinger liquid behavior depends crucially on the charge conjugation symmetry, where the Hamiltonian 
remains the same after electrons and holes are exchanged about the Fermi surface. When the valence band is 
intrinsically involved near resonance in the RRS process, such electron-hole conjugation symmetry is totally broken, 
because the filled valence band is effectively "overlapped" with the conduction band at Fermi surface. In other 
words, an electron below the conduction band Fermi surface now effectively has a new channel, not restricted by the 
small ID phase space, to be excited above the conduction Fermi surface through the two step resonant interband 
transition, through the valence band virtual transition. An estimated resonance condition for this apparent breakdown 
of LL behavior in ID RRS spectroscopy can therefore be obtained by \E^\ < qvp/Ep = 2{q/kF), which is 0.2 for 
q ~ QAkp and is consistent with our numerical result shown in Fig. 7. We therefore physically explain the failure of 
the theoretical attempt of using LL theory to study the ID RRS experiments near resonance The qualitative 

similarity of the experimental RRS results for one, two and three dimensional systems confirms our theory, which 
is based on the conventional Fermi liquid model. We note, however, that the Luttinger liquid description of ID 
systems has recently been theoretically modified in an attempt to understand the observed RRS spectra, but the 
applicable theory is quite subtle and beyond the scope of this paper. 



IV. SUMMARY 



In summary, it may be important to emphasize that the striking phenomenological similarity in the experimentally 
observed RRS spectra in one, two, and three dimensional systems is a strong indication that generic interband 
resonance physics as studied here (within a resonant RPA scheme) is playing a fundamental role in producing the 
low energy "SPE" feature in the polarized RRS spectra, which cannot be explained by the standard (nonresonant) 
theory or any other non-generic (system-dependent) theories. Our theory can also be applied to the depolarized RRS 
experiments, where both single particle and spin density excitations are important, but the exchange energy should 
be included properly to separate these two excitations which are degenerate in the regular RPA calculation. 

Once exchange correlation effects are invoked to distinguish the SPE and the SDE (with the SDE lying below the 
SPE by the exchange energy), our resonant theory can account for the observed two-peak structure in the resonant 
depolarized RRS experiments in a way very similar to the theory developed herein for the SPE and CDE in the 
polarized RRS experiments. To summarize our results, we have developed a theory for resonant Raman scattering 
spectroscopy in one, two and three dimensional semiconductor structures by considering the full two step resonance 
process involved in the scattering of external photons. Wc find that at resonance the RRS spectra have considerable 
weight at the SPE energy with the SPE weight decreasing off resonance. There is no qualitative difference in the 
RRS spectra between the systems of different dimensions. Our results are in qualitative agreement with experimental 
findings and provide a generic theoretical explanation for an ubiquitous puzzle which dates back more than twenty 
five years. As a concluding note we point out that it may be somewhat misleading to call the additional feature in the 
RRS spectra an " anomalous" SPE mode as has routinely been done in the literature — a pure SPE mode arises from 
the imaginary part of the irreducible polarizability function, as given within RPA by Eq. (8), whereas the anomalous 
additional RRS feature arises primarily from the presence of the 11^^^ term (Eq. (17)) in our resonant RPA theory 
(Eqs. (16)-(19)) which is (related to, but) quite different from the irreducible polarizability, Xlg-^^ (Eq. (8)) by virtue 
of the nontrivial nature of the resonance function A{p, q). Finally, we mention that a very recent experimental report 
has appeared in the literature |^ specifically verifying the essential features of our theory . However a complete 
quantitative understanding of experimental results may very well require inclusion of additional effects (e.g. excitonic 
corrections, many body effects) beyond the scope of our work. 

We thank A. J. Millis for critical discussion. The work is supported by the US-ARO, the US-ONR and DARPA. 
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(a) (c) 

FIG. 1. (a) Schematic representation of RRS in the direct gap two band model of electron doped GaAs nanostructure. 
and ujf are the initial and final frequencies of the external photons. Steps 1, 2 and 3 are described in the text (RRS involves 
steps 1 and 2 only) . (b) and (c) are the Feynman diagrams of the electron-photon scattering process described by p • A and 
A • A terms respectively in the interacting Hamiltonian, Hi. Solid and wavy lines represent the electron and photon Green's 
functions respectively. 




FIG. 2. Diagrammatic representation of the conduction band irreducible response function, Do ^ (q, w) and reducible 
response function, n(q,a'), in the standard random phase approximation. V'c(q) is the Coulomb interaction. 
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FIG. 3. Typical momentum- energy dispersion of the single particle excitation continuum (shaded region) and the collective 
charge density excitations (plasmons) of one, two, and three dimensional electron systems (calculated within RPA). 
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FIG. 4. Dynamical structure factor obtained by the standard (nonresonant) RPA calculation at zero temperature for (a) 

one, (b) two, and (c) three dimensional electron systems ignoring valence band effects. Solid lines are the calculated spectra 
in the long wavelength limit (small momentum transfer, |q| = O.lkp) and dashed lines are for the large momentum transfer 
(|q| = O.Tkp in (a) and (b), and |q| = OAkp in (c)) calculations. The electron densities used in the calculation is 6.5 x 10^ cm"'^, 
3.2 X 10^^ cm~^, and 1.8 x 10^^ cm~^ for one, two, and three dimensional systems respectively. Finite impurity scattering 
(7 = 1Q~^ Ef) has been included to broaden the peaks. The ripple in (b) and (c) is of numerical origin. 






FIG. 5. Diagrammatic representation of the resonant Raman scattering response function including the valence band 
electrons in the RPA calculation. Different kinds of irreducible response functions are defined and explained in Eqs. (17)- (19) 
and the matrix element, yl(k, q), is defined in Eq. (15). 
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FIG. 6. Dynamical structure factor in the resonant RPA calculation for (a) one, (b) two, and (c) three dimensional electron 
systems incorporating the valence band electrons. Wo choose the resonance broadening factor. A, to bo 0.02Ef, and finite 
impurity scattering factor, 7, to be 0.05 Ep in order to have agreement with the experimental data (the impurity broadening 
is still rather small since f/Ep = 1/20). Other system parameters are the same as used in Fig. 4. 
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FIG. 7. Ratio of the resonant Raman scattering spectral weight (CDE to SPE) as a function of the resonance energy, 
in one, two, and three dimensional systems. Off-resonance, \Eui\ > 0.2, CDE always dominates SPE in the spectra, but near 
resonance, \Euj\ < 0.2, the SPE weight could even be stronger than the CDE weight. All system parameters are the same as 
used in Fig. 6. 
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